Correlations of spin currents through a quantum dot induced by the Kondo effect 
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We study correlations of spin currents flowing through a Coulomb blockaded quantum dot. While 
vanishing for elastic co-tunneling, these correlations develop as the quantum dot enters the Kondo 
regime. They are a manifestation of Kondo physics in quantum dots. We demonstrate that the spin 
current correlator is non-perturbative in the Kondo coupling. 
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Measurements of current correlations in electrical con- 
ductors give valuable information about the microscopic 
interactions. An important example are shot noise mea- 
surements on conductors in the fractional quantum Hall 
regime pj . They have demonstrated that the excitations 
of those interacting systems carry fractions of the ele- 
mentary charge. Another well-studied interaction phe- 
nomenon in nanostructures is the Kondo effect 0]. It 
is usually observed as an increase of the conductance 
through a Coulomb blockaded quantum dot (QD) upon 
lowering its temperature Its microscopic origin are 
ground-state correlations between an unpaired spin on 
the QD and the spins of electrons in the adjacent leads. 
This suggests that a natural way to reveal Kondo corre- 
lations should be through spin-dependent quantities. In 
this letter we show that indeed Kondo physics manifests 
itself in correlations of spin currents through a QD. The 
measurement of spin currents is one of the goals of the 
rapidly developing field of spintronics Q. The applica- 
tion discussed here requires a measurement method that 
does not disrupt the Kondo effect. In particular it must 
not spin-polarize the leads. One possibility based on the 
coupling of moving spins to an electric field has been put 
forward in 

We study the correlator 



C 



Tl 



dt{i^o)h{t))-(h){h) 



(i) 



between the currents I-\ and Ii of the numbers of spin- 
up and spin-down electrons flowing through a QD. We 
assume that the applied bias voltage V and the tempera- 
ture T are well below level spacing e<2 and charging energy 
U of the dot, eV, kT < e d , U (we set h = k = 1). The 
system is then well described by the Anderson single- level 
impurity model 

H A = + H D + Un d ^n dl + vp{d\c kBa + h.c.) 
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k labels the electrons' momentum, a <E their spin 

and [3 G {L, R} the lead (left or right of the QD), n<j CT = 



d} T d a . fj,a is the electrochemical potential of lead (3. We 
have fiL — M-R = e ^ an d we choose the zero of energy at 
the Fermi level. If the dot is weakly coupled to the leads, 
<C U — e^ed (r^ = 2iri>\vp\ 2 and v is the density of 
states in the leads), and the QD is operated in a Coulomb 
blockade valley, charge fluctuations on the dot level d 
are strongly suppressed. If additionally d is occupied by 
a single electron, the spin of that electron is the only 
relevant degree of freedom of the QD. One obtains the 
corresponding low-energy Hamiltonian by a Schrieffer- 
Wolff transformation |6(, 
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are the amplitudes of scattering processes of lead 
electrons that involve the spin S of the QD, whereas 
accounts for regular potential scattering. s a — <r a /2, 
where a a are Pauli matrices. We assume that either td <C 
U - e d or t d > U - e d , such that \J^\ = \J^\. The rel- 
ative sign depends on whether transport occurs by emp- 
tying the level d, that is e d <C U - e d ( J (0) = J (0) ), or by 
doubly occupying it, that is e d 3> U — ed (J(°> — — J'- -'). 

Before giving the details of the calculation of C-j-j 
we motivate its results. For this we focus on the case 
J(°) = -J(°). In the co-tunneling regime, when only 
processes to lowest order in J are relevant, transport oc- 
curs then through virtual states with a doubly occupied 
level d. These virtual states decay into states with one 
electron of either spin on the QD, contributing to J| or I± 
with equal probabilities. To lowest order in J there are 
therefore no correlations between 7j and 2j_, C|j_ = 0. It 
is important to note, that if S — a co-tunneling docs not 
allow for the transfer of an electron with spin a without 
flipping S because of the Pauli principle. This changes 
at low temperatures, in the Kondo regime. Higher order 
tunneling processes can then transfer electrons with ei- 
ther spin without flipping S. Being equally likely, these 
"non-spin- flip" processes do not contribute to C^. The 
spin current produced by processes that flip S, however, 
is correlated with the state of S. Whenever the spin on 
the dot is "up" , it can flip down and make a contribu- 
tion to If. If S is "down" , a spin-flip process produces a 
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FIG. 1: Illustration of the mechanism that introduces spin- 
current correlations in the Kondo regime. Only the relevant 
spin-flip scattering processes are depicted. The QD has dif- 
ferent spin at consecutive scattering times t and t + At. Con- 
sequently, the outgoing electrons carry alternating spin. 



tended to non-equilibrium situations as considered here. 
Following [9| we eliminate virtual transitions into high 
energy electron states in the leading logarithm approx- 
imation down to a scale £1 > eV. This results in an 
effective Hamiltonian H that has the form of Eq. J3J), 
but renormalized coupling constants J and J instead of 
and and reduced bandwidth A (instead of A^) 
with e(A) = Q. Spin fluctuations enhance the spin- flip 
amplitude to J = l/2i/ln(i7/Tj<-), where Tk is the Kondo 
temperature, while J — Our limit vJ <C 1 implies 

eV » Tk- We additionally assume that the tempera- 
ture is smaller but of the same order of magnitude as the 
voltage, eV > T. 

The central object in our approach is the generating 
function of spin-current correlators 



spin-down current. Since in the absence of external spin- 
relaxation processes every flipping of S to "up" has to 
be followed by a spin- flip to "down" , both processes oc- 
cur alternatingly and equally often, as illustrated in Fig. 
^ This introduces correlations between if and if and 
we expect Cff to be non-vanishing in the Kondo regime. 
Formally, the Pauli blocking in the co-tunneling regime 
is described by Hamiltonian through an interference 
of the amplitudes and . It is lifted in the Kondo 
regime that is described by an effective Hamiltonian of 
the form © with > 

We study the limit of a weak Kondo effect with small 
effective low-energy couplings vJ <C 1. One might ex- 
pect that perturbation theory in J is then adequate to 
calculate Cff. This perturbation expansion is, however, 
plagued by divergences due to the fact that the corre- 
lation functions of S do not decay in time to low or- 
der in J. This indicates that long-time spin correlations 
have to be treated non-perturbatively. For this we use 
a method that is close in spirit to the master equation 
approach to current correlations in sequential tunneling 
by Bagrets and Nazarov [3. Since we are interested in 
the co-tunneling and the Kondo regime we cannot di- 
rectly apply their method and we derive a variant of it. 
Our approach is quantum mechanical at short time scales 
while non-perturbative and classical at times that exceed 
the decoherence time of S. 

In the limit of interest v J <C 1, we can apply the per- 
turbative renormalization group [8j to capture the short- 
time spin dynamics. In 9] this approach has been ex- 



Z{\) = Tre-*"" 1 " e~ im e l X ^ N " e lHt . (4) 

p( m > is the initial density matrix of the conductor and 
Z generates moments of the number AN^ a of spin a 
electrons that is transferred into lead f3 during time t, 
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Z as defined in Eq. is a Keldysh partition function. 
Its two time development operators can be implemented 
by operators on the two branches of the Keldysh time- 
contour 0. The exponentials exp{±i\p cr N l3 ' J } can be 
absorbed into phases for the tunneling terms in H. By 
inserting complete sets of Fermion coherent states [llj we 
then derive a mixed representation |l2j | of Z as a path 
integral over a time-ordered spin operator expression, 



Z{\) =Tr s T c 
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The Fermion fields c carry Keldysh, lead, spin, momen- 
tum, and frequency indices a, (3, a, k, and lu. Q, a matrix 
in this space, is the electron Green function correspond- 
ing to Hl- J x describes tunneling processes and con- 
tains spin operators. T c denotes operator ordering along 
the Keldysh time-contour and relative to the initial spin 
density matrix pi . The trace Trs over spin operators 
is taken. We have 



with the third Pauli matrix t z and spin operators S"' uu = J dte^ 1 ^^^ )*5" (t) on the Keldysh branch a. Integrating 
over the Fermions in Eq. © we obtain 

Z(X) =Tr s T cP P e Tli ^ 1+GJ ^ Z\ J=J=0 = Trs T c pP e T ^ gjX -^ jXgjX+ °^ = Tr s T c pp } e~ c \ (8) 
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The 0{J) term of C x does not contain spin operators, while the terms of 0(J 3 ) are negligible for vJ -C 1. To save 
space we evaluate here only the 0(J 2 ) term C x that derives from the spin off-diagonal matrix elements of J x , 

£0= E J^T^^t-A^M.^ya I dtd£GZ^-t)G%{t-l!)^{t)^), (9) 
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where G^ a ,(t) = / (dw/27r) e luJt Qa/3tku,a> p\ku and 
5" = S" ±iSy. The Green functions in Eq. (0 summed 
over fc, A:' decay exponentially over the time tt = 1 /T. At 
time scales longer than tt £q is therefore local in time. 
It moreover couples spin-flip processes occurring at the 
same time on different branches of the Keldysh contour. 
This leads to classical behavior at long times and it al- 
lows for a description along the lines of Q. One could 
obtain the effective theory on the scale tt by integrating 
out fast spin fluctuations in a path integral for the orig- 
inal model Eq. (JSJ 13]. Equivalently we integrated out 
high frequency spin fluctuations in the Hamiltonian for- 
malism. Due to the form of the Green functions in the 
effective model H with reduced bandwidth, C x is then 
strongly suppressed for frequencies lu > Q and these fluc- 
tuations contribute to Z only negligibly. The corrections 
to the scaling logarithms due to spin fluctuations in the 
frequency range T < lu < f2 are for our choice of pa- 
rameters f2 > eV > T of O(l) and to leading order in 



the logarithms negligible as well. To this accuracy C x 
therefore equals the corresponding piece of the effective 
theory at the scale Tt- The contribution due to the diag- 
onal elements of J x has an analogous structure and we 
conclude that C x is local on the time scale tt ■ In Eq. JSJ 
we have expressed Z as the trace over a time dependent 
density matrix pg(t) = T c p^ exp(— C x ). Because of the 
locality of C x , p^ obeys an ordinary differential equation. 
The off-diagonal entries of pg decay exponentially under 
evolution with that equation. This shows that the spin 
dynamics is indeed classical on the time scale Tt- It is 
therefore sufficient to study the evolution of the diagonal 
elements p x and p x of p s , that we collect into a vector 
p x = (p x ,p x ). p x obeys 

d tP x = -L x p x , (10) 
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For conciseness we assume from now on that the QD is 
symmetrically coupled to the leads, Jpp' = J, Jpp = J. 
a/3j3'\J\ 2 are the usual rates for tunneling between two 
leads. Ts — Yla/3' a f3t3'\J\ 2 + 7S is the relaxation rate of 
S. The first term accounts for spin- flips by conduction 
electrons. The second term 75 has been introduced phe- 
nomenologically. It accounts for spin-flip processes that 
are not described by our model. We assume 7s <C T, 
such that this additional spin-decoherence does not af- 
fect Kondo correlations. Eqs. (jlOJI with (fill) differ from 
their counterparts for sequential tunneling |jj mainly in 
how the electron state occupation numbers enter. Kondo 
correlations are accounted for by the renormalization of 
J. The off-diagonal entries of L x describe spin-flips of S. 
Because electrons can also be transferred without flipping 



S, the diagonal elements of L x contain counting factors 
C x as well. 

We integrate Eq. IjlUfl to obtain Z. At long times 
t it simplifies to the exponential of the smallest eigen- 
value of —tL x . C-\i can then be obtained using Eq. JSJ. 
We introduce transmission probabilities r = (2-kvJ) 2 , 
f = (2-kvJ) 2 and the spin-off-diagonal Fano-factor — 
On/I (I is the mean current through the QD). Although 
not within the limits of applicability of our theory, it is 
instructive to first take the zero temperature limit 
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Eq. 113J1 displays most clearly the main finding of this 
letter: While vanishes in the absence of Kondo cor- 
relations, when potential and spin-flip scattering pro- 
cesses have equal probabilities r = f, it grows non-zero 
in the Kondo regime. We need to convince ourselves 
that this carries over to finite temperature T ~ eV. 
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The first correction to Eq. I|13|) due to temperature, 
Af' 1 ' = [(f - t)/(3t + f)](T/eV) (7s = 0) is reassur- 
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FIG. 2: Voltage dependence of the Fano factor Fji for T — 
eV/30, e d = 1 meV, and T K = 10" 7 e d (solid line: 7s = 0, 
dashed line: 7s = 10 J s _1 ). Ffi develops in parallel with 
Kondo correlations upon lowering the voltage/temperature. 

is shown in Fig. [21 In this plot both temperature and 
voltage are varied, at constant ratio eV/T > 1. Although 
the Fano factor at finite temperature does not reach 
its theoretical zero temperature maximum 1/6, its qual- 
itative behaviour is robust: it develops in parallel with 
Kondo correlations. The effect is cut-off at small volt- 
ages eV <C 7s/t by external spin-flip processes. S is 
then flipped randomly in between electron transfers and 
the mechanism illustrated in Fig. ^ is inoperative. For 
Fig. |2| we have assumed typical experimental parameters 
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e c i = 1 meV and 75 = 10 5 s _1 . Even longer spin relax- 
ation times have been observed [l4| . 

We come back now to the difficulties encountered with 
perturbation theory in J. The zero temperature limit 
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cannot be accessed in an expansion around J = 0. It lies 
outside its radius of convergence. Perturbation theory in 
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by a straightforward extension of our method to slowly 
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The dispersion of Cj| on the scale Tg, that is of the order 
of the mean current I/e at small 7s, is a consequence of 
the time-ordering of spin-flips described above. 

In conclusion, we have studied correlations of spin cur- 
rents through a QD. We have found that these correla- 
tions can be induced by Kondo fluctuations on the QD. 
They are a new manifestation of the Kondo effect in QDs. 
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